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In view of the wide success of molecular quasi-chemical theory of liquids, this paper develops the 
soft-cutoff version of that theory. This development has important practical consequences in the 
common cases that the packing contribution dominates the solvation free energy of realistically- 
modeled molecules because treatment of hard-core interactions usually requires special purpose 
simulation methods. In contrast, treatment of smooth repulsive interactions is typically straightfor- 
ward on the basis of widely available software. This development also shows how fluids composed 
of molecules with smooth repulsive interactions can be treated analogously to the molecular-field 
theory of the hard-sphere fluid. In the treatment of liquid water, quasi-chemical theory with soft- 
cutoff conditioning doesn't change the fundamental convergence characteristics of the theory using 
hard-cutoff conditioning. In fact, hard cutoffs are found here to work better than softer ones. 



I. INTRODUCTION 

Recent developments of quasi-chemical theory consti- 
tute a fresh attack on the theory of solutions f rom a 
molecular scale, a point that has been noted before 
Quasi-chemical theory provides a natural organization 
for calculations of ion-water chemical int eract ions in the 
treatment of ion hydration free energies ! 4 I 5 I 6 I A remark- 
ably natural quasi-chemical theory for the statistical 
thermodynamics of liquid water^ can be implemented on 
the basis of the massive data sets of molecular simulations 
and provides a realistic direct evaluation of the entropy of 
that liquid. For non-polar solutes dissolved in liquid wa- 
ter, quasi-chemical theory avoids physical assumptions of 
the van der Waals (perturbative) type that are otherwise 
common, and thus resol ves questions centering on refer- 
ence system assumptions!^ As originally developed with 
a sharp cutoff to spatially define the inner shell, quasi- 
chemical theory implemented as a molcculc-ficld theory 
for the hard-sphere fluid provides an equation of state as 
accurate as the most accurate previous theory.^ 

This quasi-chemical theory can be expressed entirely as 
statistical modeling of the distribution P(e) of energies e 
binding a distinguished molecule to the medium. A goal 
is to discriminate effects of close, or inner-shell, neigh- 
bors of that molecule from the effects of the more distant 
material. This is done by formulating a statistical con- 
dition to identify inner-shell molecules. We will express 
our condition with \, a function indicating that the inner 
shell is empty when x = !• X nas been assumed to be ei- 
ther or 1, sharply changing between those possibilities. 
We then direct attention to the conditional probability 
P(e\\ = 1) in order to define an outer-shell contribution, 
and then further use chemical concepts to evaluate the 
remainder associated with occupancy of the inner-shell. 



This quasi-chemical theory with sharp- or hard-cutoff 
conditioning works satisfactorily both for the free en- 
ergy of liquid watei^ and for the hydration free energy of 
CF4(aq)Pln the latter case, the theory generates a re- 
sult similar to a van der Waals theory without requiring 
those assumptions in advance. A prominent contribu- 
tion to that result is a hard-core or packing contribution 
associated with the sharp \ = 1 condition. On an in- 
tuitive basis, it is widely asked whether a more natural 
packing contribution might be based upon smooth repul- 
sive interactions. In the contrasting case of liquid water, 
more aggressive conditioning is required, and the pack- 
ing contribution clearly does not have the logical status 
of a contribution of repulsive-force reference interaction 
considered to be similar to actual intermolecular inter- 
actions. The present work is designed to clarify these 
contrasting points. 

We anticipate the results below by noting here that 
the theory for the soft-cutoff conditioning can be devel- 
oped in straightforward analogy with that for hard-cutoff 
conditioning. This development has important practical 
consequences because a packing contribution frequently 
dominates the solvation free energy of realistically- 
modeled molecules but treatment of hard-core interac- 
tions associated with a hard cutoff usually requires spe- 
cial purpose calculations. Examples of the special con- 
siderations that apply to simulation of liquids of hard- 
core-model polyatomic molecules can be found in the 
references! 11 * 12 ^ In contrast, treatments of smooth repul- 
sive interactions are typically straightforward on the ba- 
sis of widely available software. The present develop- 
ment implies also that fluids composed of molecules with 
smooth repulsive interactions can be treated by analogy 
with the molecular-field theory of the hard-sphere fluid P^ 

In a case such as CF^aq) where only minimal condi- 



2 



tioning is required,^ we propose below a soft cutoff that 
should reduce the variance of estimates of the required 
contributions. In the case of liquid water where aggres- 
sive conditioning is required,^ we can systematically vary 
the softness of the cutoff defining the inner shell, and at- 
tempt to learn what works best. 



II. THEORY 

Quasi-chemical theory with sharp- or hard-cutoff con- 
ditioning, which we aim to generalize, can be expressed 
as 

/3/i cx = - In «x)) +ln ( X )+ln / e? e P(e\ X = l)de , (la) 



/? M cx «-ln((x» + ln( X > 

+ f3(e\ X =l) + y <& 2 |X = 1) • (lb) 

The notation ((. . .)) indicates the average over the ther- 
mal motion of the system and the solute with no inter- 
action between them, thus the subscript 0. x is an indi- 
cator function and x = 1 defines our condition. As an 
example, for the application to liquid water x = 1 when 
there are no O-atoms of bath molecules within a radius A 
from the O-atom of a distinguished water molecule. Then 
X = if any solvent molecule is present in that inner- 
shell. The conditional averages indicate that statistics 
are collected only when x — 1- Significant conditioning 
produces slightly sub-gaussian binding energy probabil- 
ity densities P(e\x — !)■ Slight sub-gaussian behavior 
leads to errors of 1 kcal/mol in the excess chemical poten- 
tials calculated with Eq. (flbl) IZ113J Nevertheless, Eq. (flbl 



or Eq. (la I after evaluating the logarithm. For a sharp 



adopts the gaussian approximation. 

The generalization that we seek can be economically 
identified on the basis of a few primitive relations. The 
first is the rule of averages^ 

(F) = < t"-Trm )) ° - ^ ((^ AU F))o ■ ( 2 ) 



Then, 
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is an example of specific utility for our argument. For an 
indicator function x we have in addition 

(G X ) = (G\x = 1) (X) , (4) 
and (x) is a probability. Collecting these relations yields 



(5) 



cutoff case, this is obtained directly from the stan- 
dard thermodynamic difference formula.^^ Beyond tra- 
ditional statistical thermodynamics, however, this for- 
mula only depends on the fact that x is an indicator 
function. 

We obtain soft-cutoff conditioning by designing an in- 
dicator function which categorizes a solvent molecule at 
a particular configuration to be either in or out of the in- 
ner shell according to a statistical rule. For example, let 
the designated or solute molecule be at a fixed position, 
and consider the i th solvent molecule again at a definite 
configuration denoted by 72,. We take < b(lZi) < 1 
to be the probability that the i th solvent molecule is 
in, independent of any other characteristic. [For the 
cases considered above 6(72^) only depends on the three- 
dimensional position of the O (oxygen) atom of the i th 
water molecule.] For each configuration, we could then 
sample b(JZi) for each solvent molecule, and thereby iden- 
tify each solvent molecule as either in or out. Exhaustive 
in-out sampling at a fixed configuration produces 
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as a sampling-averaged indicator function. The formula 
specifically indicating a soft-cutoff is then 



<e^| X =l) 



( (x) 



V((X» 



(7) 



The hard-cutoff theory is recovered if step functions are 
chosen for the b(lZi). 

It is interesting to note that Eq. Q is specialization of 



(e^ u W) 
{{W)) 



(8) 



to the case that W — X i s an indicator function. 
BennettP^ has studied such forms without the restric- 
tion that W be an indicator function, and found W oc 
(l+ze' 3A!7 ) _1 to be statistically optimal with a particular 
value of z. In contrast, our goal is the dis-entanglement 
of the effects of longer-ranged intermolecular interactions 
from packing and chemical contributions, so we do not 
follow Bennett's results here. A potential advantage of 
Eq. is that (x) is an occupancy probability and there- 
fore can be modeled on a pattern of chemical equilibria, 
as we elaborate below. 

Nevertheless, if minimal conditioning is all that is re- 
quired, then a straightforward suggestion for choice of the 
indicator function is x — exp (— 0uq) for a smooth repul- 
sive molecular-pair interaction, uq, that matches AU at 
short range. In view of Eq. (J8|, this choice would re- 
duce the variance of the numerator of Eq. ([T]) . Then the 
denominator is ((x))o = ex P ( — /^o*)- 
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Chemical interpretation 

A general derivation of quasi-chemical theory 15 can be 
based upon an inclusion-exclusion development initiated 
from the tautology 



i=np'(^) + ( 1 - 6 ^))] 



(9) 



When included within the averaging brackets that yield 
the chemical potential, Eq. [9] leads to the general 
development)^ 



l + J2K n p n ]x((e-^x)) ■ (10) 



Since 
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Eq. (10) implies 



<x> = 



1 + £ K n p" 

n>l 



(11) 
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The K n 's involve the cutoff function in a natural way™ 
Together with 



p(n) 




n w 

3=1 



' k=N 

n (i-w) 

A-n+l 



,(13) 



K n p n is the ratio of the probability that the distinguished 
species has n ligands to the probability that it has zero 
(0) ligands 



p(0) 



K n p n 



(14) 



This is naturally interpreted on the basis of the chemical 
equilibrium 



W* 



nW ^ w*w„ 



(15) 



where W* denotes a water molecule with zero (0) ligands 
and W*W„ denotes that water molecule with n ligands 
occupying the inner shell. Given Eq. (12 1, the probabil- 



ity that the distinguished molecule has zero (0) ligands 
is (x). This establishes the thermodynamic significance 
of that normalization factor. These interpretations are 
straightforward for the sharp-cutoff case. If the proba- 
bilities of Eq. (l4"| are given by a Poisson distribution, 
then K n — v n /nl with v the volume of the stencil. For a 
specified function b{1Zi), p{n) and K n can be observed in 
conventional simulations where the full physical solute- 
solvent interactions are expressed. 



These analogies go further. For example, consider the 
packing contribution ((x)) J^and apply Eq. (10 1 to the 
case AU a — 0. This yields 



or 



i = i + £ k nP n x « X )) 



«X»o 



1 + Y, K np H 
n>l 



(16) 
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with 



P{n) 
p(0) 



K n p" 



(18) 



Here the tildes indicate quantities observed without the 
solute physically present, but with a soft stencil defined 
by the indicator function \. Specifically 




n 



3 = 1 



k=N 



n (i-w) 



.k—n+l 



(19) 



i.e., Eq. (13 1 for the present uncoupled case. Fig. [T] 
gives example results of p (n) for liquid water, and shows 
that soft-cutoff conditioning produces integer-referenced 
probabilities. 
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FIG. 1: For TIP3P water at T=300K and p=latm: p(ri) 
corresponding to the soft-cutoff function of Eq. ( 23 1 with A = 
3.8A. The smooth curves plotted are two moment maximum 
entropy models obtained on the basis of the Poisson default 
modelP 



If we take \ = ex P ( — P u o) as suggested above, Eq. ( 17 1 
j;ives the new and surprising result 



= 1 + E^ 



p 



(20) 



n>l 
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for the free energies of model liquids composed of 
molecules with soft repulsive interactions. Eq. (20) is, 



however, directly analogous to the quasi-chemical the- 
ory for hard-core fluids. 111 ^ The utility of Eq. |20) is 
due, in part, to the fact the low-density limiting values 

fim p ^o K n = Kn(T) can be evaluated once-and-for-all, 
and they provide a re mark ably natural initial model for 
that equation of stateA^ Collecting preceding results, 
Eq. can be expressed in the appealing form 



^"" = (e" Atr | X =l) 



'1+ E K n p n 

n>l 

1 + E K n p" 

n>l 



(21) 



The developments of this subsection are expected to 
have the practical consequence that the whole of theory 
of Ref . 8 can be consistently implemented on the molecu- 
lar dynamics simulations as contrasted with Monte Carlo 
calculations. 



Though results such as Eq. ( 20 1 give a compelling phys 



ical organization to the statistical problem of evaluating 
we emphasize that this approach is thoroughly tau- 
tological. For example, the identification Eq. (18 1 trans- 
forms Eq. pOj) to 



^p(O) 

which is transparently correct. 



1 

W) 



(22) 



III. COMPUTATIONAL IMPLEMENTATION 

All calculations were performed with the MUSIC 
softwartpEH which is a set of FORTRAN 90 modules 
for MC and MD simulations. A cubic simulation cell of 
edge- length 30 A was used in all NVT Monte Carlo simu- 
lations. Water was modeled with the TIP3P forcefieldP^ 
Based on the density of TIP3P water at 300K and 1 
atmJiSl 888 water molecules were used in the simula- 
tions. All Lennard- Jones and electrostatic interactions 
were cutoff at 9.0 A. The electrostatic interactions were 
described by the shifted force method of Fennell and 
GezelterP 

From a simulation record, the excess free energy can 
be directly evaluated by applying Eq. Specifi- 
cally for the soft cutoff, averages such as (x), (ex) = 
(e | x = 1) (x)t an d so on, can be evaluated straightfor- 
wardly. These were obtained for the smooth cutoff func- 
tion 



1 - b(r) = 



1 



l + e -(r-\)/R ' 



(23) 



with R a softness parameter. The results of applying 
soft-cutoff for water at 300K are shown in FIG. [21 The 
results there are all qualitatively similar in overshooting 
the numerically exact value to which they must return 



with A — > oo. That overshoot is an indication that the 
gaussian model places slightly too much weight in the 
high (unfavorable) e tailpESj.e., the actual distribution 
is sub-gaussian. 




FIG. 2: The free energy discrepancy on the basis of Eq. ( lb I 



for TIP3P water at 300K and normal density utilizing differ- 
ent smooth cutoff functions of the form Eq. (23 1. 



With the present generalized theory in hand, we in- 
vestigated the following idea: Consider a family of cutoff 
functions for which the softness can be systematically 
varied. An example is 



1 - b m (r) 



-(A/r) 



(24) 



with a parameter m > 3. With increasing m this cut- 
off becomes increasingly sharp. A natural idea then is 
to choose m to minimize the variance of the conditional 
distribution P(e\x — 1) appearing in Eq. (la I. The goal 
is to narrow the distribution so that deviations from the 
gaussian model for /x ex (Eq. ( lb I ) are less important. We 



found that the variance decreases with increasing m (in- 
creasing sharpness) towards an m ~ oo limiting value, 
the hard-cutoff case as a practical matter. Hard-cutoff 
conditioning is always preferred for such cases. The con- 
clusion is that for systems and calculations such as these, 
the results are always improved the more that the con- 
ditioning eliminates statistical modeling of short-ranged 
interactions, and hard cutoffs do that best. Those in- 
teractions include both associative and repulsive interac- 
tions that encompass both chemical and packing contri- 
butions. This conclusion is clearly consistent with the 
results shown in FIG. [2j 



IV. CONCLUSIONS 

Quasi-chemical theory with a soft-cutoff can be de- 
veloped in straightforward analogy with the hard-cutoff 
theory. This development has important practical conse- 
quences in the common cases that the packing contribu- 
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tion dominates the solvation free energy of realistically- 
modeled molecules because treatment of hard-core inter- 
actions usually require special purpose simulation meth- 
ods. In contrast, treatment of smooth repulsive interac- 
tions is typically straightforward on the basis of widely 
available software. This development also shows how flu- 
ids composed of molecules with smooth repulsive inter- 
actions can be treated analogously to the molecular-field 
theory of the hard-sphere fluid™ 

In the treatment of liquid water, soft-cutoff condition- 
ing doesn't change the fundamental convergence charac- 



teristics of the hard-cutoff approach. Hard cutoffs are 
found to work better than softer ones in that case. 
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